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We explore the cosmological implications of an ultra-light pseudo-Nambu- 
old 

10 18 GeV and explicit breaking scale comparable to MSW neutrino masses, 



q Goldstone boson. With global spontaneous symmetry breaking scale / 

h 

M ~ 1(T 3 eV, such a field, which acquires a mass ~ M 2 / f ~ H , would 
have become dynamical at recent epochs and currently dominate the energy 
density of the universe. The field acts as an effective cosmological constant 
for several expansion times and then relaxes into a condensate of coherent 
non-relativistic bosons. Such a model can reconcile dynamical estimates of 
the density parameter, Q m ~ 0.2, with a spatially flat universe, and can 
yield an expansion age H^tQ ~ 1 while remaining consistent with limits from 
gravitational lens statistics. 
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Recently, a cosmological model with substantial vacuum energy — a relic cosmological 
constant A — has come into vogue for several reasons. First, dynamical estimates of the 
mass density on the scales of galaxy clusters, the largest gravitationally bound systems, 
suggest that £l m = 0.2 zb 0.1 for the matter (m) which clusters gravitationally (where the 
density parameter Q is the ratio of the mean mass density of the universe to the critical 
Einstein-de Sitter density, tt(t) = 8?rGp/3i? 2 ) [1]. H owever, if a sufficiently long epoch 
of inflation took place during the early universe, the present spatial curvature should be 
negligibly small, fl tot = 1. A form of dark, homogeneously distributed energy density with 
Clh = 1 — Cl m , such as a cosmological constant, is one way to resolve the discrepancy between 
Vt m and fltot- 

The second motivation for the revival of the cosmological constant is the 'age crisis' for 
spatially flat fl m = 1 models. Current estimates of the Hubble expansion parameter from a 
variety of methods, most recently Cepheid variable stars in the Virgo cluster [2], are (with 
some notable exceptions) converging to relatively high values, Hq ~ 80 ± 15 km/sec/Mpc [3], 
while estimates of the age of the universe from globular clusters are holding at t gc ~ 13 — 15 
Gyr or more [4]. Thus, the 'expansion age' HqIq = 1.14(i?o/80km/sec/Mpc)(to/14Gyr) 
is uncomfortably high compared to that for the standard Einstein-de Sitter model with 
fl m = 1, for which HqIq = 2/3. On the other hand, for models with a cosmological constant, 
i?o^o can be significantly larger: for example, for = A/3i?Q = 0.8 = 1 — fl m , one finds 
HqIq = 1.076. Third, cosmological constant-dominated models for large-scale structure 
formation with cold dark matter (CDM) and a nearly scale-invariant spectrum of primordial 
density perturbations (as predicted by inflation) provide a better fit to the observed power 
spectrum of galaxy clustering than does the 'standard' fl m = 1 CDM model [5]. 

While they provide a number of theoretical benefits, models with a relic cosmological 
constant have problems of their own. A cosmological constant for which, e.g., fl\ ~ 1 
corresponds to a vacuum energy density p vac = A/8vrG ~ (0.003 eV) 4 . Within the context 
of quantum field theory, there is as yet no understanding of why the vacuum energy density 
arising from zero-point fluctuations is not of order the Planck scale, Mp h or at least of 
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order the supersymmetry breaking scale, Mg USY ~ TeV 4 , both many orders of magnitude 
larger. Within the context of classical field theory, there is no understanding of why the 
vacuum energy density is not of the order of the scale of one of the vacuum condensates, 
such as .\/ ( ' ; , ,. -Mj USY , -M^sm 4 e w /(4TTa) 2 ~ (175 GeV) 4 , or -/ 4 - (100 MeV) 4 . 
Thus, a vacuum density of order (0.003 eV) 4 appears to require cancellation between two 
(or more) large numbers to very high precision. In addition, it implies that we are observing 
the universe just at the special epoch when fl m is comparable to fl\, which might seem to 
beg for further explanation. 

Moreover, such models now face strong observational constraints from gravitational lens 
statistics: in a spatially flat universe with non-zero A, the lensing optical depth at moderate 
redshift is substantially larger than in the Einstein-de Sitter model with fl m = 1 [6]. In 
the Hubble Space Telescope Snapshot Survey for lensed quasars, there are only four lens 
candidates (thought to be lensed by foreground galaxies) in a sample of 502 QSOs; from this 
data, the bound fl\ ^ 0.6 — 0.8 has been inferred [7]. For fl\ = 1 — Qq < 0.7, the expansion 
age satisfies i?oto < 0.96. With a cosmological constant saturating this bound, the globular 
cluster age to > 14 Gyr implies Hq < 67 km/sec/Mpc, within the uncertainties of but below 
the central value of recent Hubble parameter determinations. 

It is conventional to assume that the fundamental vacuum energy of the universe is zero, 
owing to some as yet not understood mechanism, and that this new physical mechanism 
'commutes' with other dynamical effects that lead to sources of energy density (after all, 
there is gravitational energy density acting on cosmological scales). This is required so that, 
e.g., at earlier epochs there can temporarily exist non-zero vacuum energy which allows 
inflation to take place, but the situation in reality could be more complex. Nonetheless, if 
this simple hypothesis is the case, then the effective vacuum energy at any epoch will be 
dominated by the heaviest fields which have not yet relaxed to their vacuum state. At late 
times, these fields must be very light. This is a big asumption: the cosmological 'constant' 
may be in the process of relaxing in a self-consistent way which leaves a residual effect at 
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any scale, and we can only hope that this hypothesis approximates this possibility. 

Adopting this working hypothesis, in this Letter we explore the consequences of an 
ultra-light pseudo-Nambu-Goldstone boson (hereafter, PNGB) field which is (i) currently 
relaxing to its vacuum state and which (ii) dynamically dominates the energy density dur- 
ing the epoch in which it relaxes. PNGB models are characterized by two mass scales, a 
spontaneous and an explicit symmmetry breaking scale; we will see that the two dynami- 
cal conditions above essentially fix these two mass scales to values which are 'reasonable' 
from the viewpoint of particle physics. Since these scales can have a plausible origin in 
particle physics models, we may have an explanation for the 'coincidence' that the vacuum 
energy is dynamically important at the present epoch. Moreover, in these models, the cos- 
mological constant is evanescent, within a few expansion times converting into scalar field 
oscillations which subsequently redshift as non-relativistic matter. Thus, unlike cosmologi- 
cal constant-dominated models, the universe is not now entering a phase of exponential de 
Sitter expansion, but has rather undergone a brief hiatus of quasi-accelerated expansion. 
As a byproduct, we shall see that the gravitational lens constraints on HqIq in this model 
are slightly less severe than for cosmological constant models, allowing an expansion age as 
large as H t ~ 1.05. 

In particle physics, the best known example of a PNGB is the ordinary 7r meson (the 
longitudinal W and Z bosons are actually exact Nambu-Goldstone bosons in association 
with gauge fields). An example of a very light hypothetical PNGB is the axion, associated 
with the Peccei-Quinn symmetry introduced to solve the strong CP problem [8]. Axions arise 
when a global U(1)pq symmetry is spontaneously broken by the vacuum expectation value 
of a complex scalar at the scale f a , ($) = f a e ia ^ a ; at this scale, the axion, the angular field a 
around the infinitely degenerate minimum of the potential, is a massless Nambu-Goldstone 
boson. QCD instantons explicitly break the global symmetry at the scale f T ~ 100 MeV, 
generating the axion mass, m a ~ 0(m ff / T // a ). Since its couplings and mass are suppressed 
by inverse powers of f a , the axion is very light and very weakly interacting. Nevertheless, 
it can play an important role in astrophysics and cosmology; indeed, astrophysical and 
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cosmological arguments constrain the global symmetry breaking scale to lie in a narrow 
window around f a ~ 10 10 — 10 12 GeV. Thus, the axion mass m a ~ 10~ 5 eV(10 12 GeV// a ), 
and its Compton wavelength is macroscopic, A a ~ (/ a /10 12 GeV) cm. 

Although motivated by the strong CP problem, the axion is a particular instance of a 
more general phenomenon that includes familons, majorons, [9] and more exotic objects [10]. 
In all these models, the key ingredients are the scale of spontaneous symmetry breaking 
/ (at which the effective Lagrangian still retains the symmetry) and a scale of explicit 
symmetry breaking \i (at which the effective Lagrangian contains the explicit symmetry 
breaking term). The mass of the PNGB is then ~ yU 2 //- Ref. [11] introduced a class of 
PNGBs closely related to familons (called 'schizons'), with masses ~ m2 f er mionl f ■ Models 
in which m fermion is associated with a hypothetical neutrino mass, m v ~ 0.001 — 0.01 eV, 
and / ~ Mqut — Mpi ~ 10 15 — 10 19 GeV, were studied in ref. [12] in the context of late 
time phase transitions [13] and form the theoretical basis for the present work. In this case, 
the PNGB Compton wavelength m^ 1 is comparable to cosmological distance scales. 

From the viewpoint of quantum field theory, pseudo-Nambu-Goldstone bosons are the 
only way to have naturally ultra-low mass, spin-0 particles. In this regard, 'technically' 
natural small mass scales are those which are protected by symmetries, such that when the 
small masses are set to zero, they cannot be generated in any order of perturbation theory, 
owing to the restrictive symmetry. For generic PNGBs, when the symmetry breaking scale 
yu is set to zero, the symmetry becomes exact, and radiative corrections do not yield an 
explicit symmetry breaking term (the radiative corrections are "multiplicative" of the scale 
li in this situation). In the ultra-light PNGB models mentioned above, the small mass 

is protected by fermionic chiral symmetries (and additional discrete symmetries) and is 
therefore technically natural. That is, when certain fermion mass terms are set to zero in 
the Lagrangian, the PNGB mass goes to zero; the fermion mass terms will not be generated 
in any order of perturbation theory. 

As an example, consider the Zjy-invariant low-energy effective chiral Lagrangian for N 
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neutrinos [12], 

-i N-l 

C = -d^cp + Y, VjiYdM + (mo + ee^ +2 ^) v ]L v ]R + h.c. (1) 
z j=0 

where f(R t L) are respectively right- and left-handed projections, V{r,l) = (1 ±7 5 )z^/2. The 
term proportional to e can arise from a Yukawa coupling qvlVrQ + h.c, where the complex 
scalar field $ has a non-zero vacuum expectation value, {$} = fe l< ^^/\^2, and e = gf 
The term proportional to mo is an explicit breaking which usually comes from some deeper 
breaking in the theory In the limit m$ — >• 0, this is a familiar chiral Lagrangian, possessing 
a continuous £7(1) chiral symmetry. The £7(1) chiral symmetry is broken to a residual Zjy 
discrete symmetry: 

vj -> z/,+1 ; z/ w _i ^ z^o J 4>^> (f) + 27rjf/N . (2) 
The induced one-loop correction, with cutoff A < /, is 

A-z«* = E T^ ln (^2) . ( 3 ) 

where 

M? = m 2 + e 2 + 2m e cos (^ + ^pj , (4) 

which respects the discrete symmetry. For N = 2, the leading contribution is log divergent, 
and the induced PNGB mass is of order ~ m^ej f\ if e ~ mg ~ m„, then ~ m 2 / f. 
For A^ > 2, the sum YljMj is independent of thus, the ^-dependent term is independent 
of the cutoff A, and for N > 2 we can write the ldoop effective potential, 

M 4 

n<f>) = -Y^Mf. ( 5) 

In this case, the ^-potential is explicitly calculable, and one again finds a quasi-periodic 
potential with mass scale ~ m 2 ,/ f . 

We are thus led to study the cosmological evolution of a light scalar field <f> with effective 
Lagrangian 

C = U^d^ - M 4 [cos(0//) + 1] . (6) 
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The theory is determined by two mass scales, M, which from (1) is expected to be within 
an order of magnitude of a light fermion (neutrino) mass, and /, the global symmetry 
breaking scale. Since <f> will turn out to be extremely light, we assume that it is the only 
classical field which has not yet reached its vacuum expectation value. Thus, in accordance 
with our working hypothesis, the constant term in the PNGB potential has been chosen to 
ensure that the vacuum energy vanishes at the minimum of the (f> potential. We focus upon 
the spatially homogeneous, zero-momentum mode of the field, <f>(t) = {(j)(x,t)), where the 
brackets denote spatial averaging. We are assuming that the spatial fluctuation amplitude 
6<p(x,t) is small compared to (f)(t), as would be expected after inflation if the post-inflation 
reheat temperature Tr# < /: in this case, aside from inflation-induced quantum fluctuations 
(which correspond to isocurvature density perturbations [15]), the field will be homogeneous 
over many present Hubble volumes. Since we will be interested in the case / ~ Mpi (see 
below), this is not a significant restriction. Finally, for simplicity we assume that any 
finite-temperature corrections to the potential V(<f>) in (6) are unimportant at the epochs 
of interest (this is different from the case of axions, for which finite-temperature corrections 
do affect the axion field evolution). The scalar equation of motion is then 

4> + 3H<j) + dV(<f>)/d<f> = , (7) 

where the Hubble parameter is given by H 2 = (a/a) 2 = (87r/3Mp ; )(p TO + p^) for a spatially 
flat universe, fl m + = 1, a(t) is the cosmic scale factor, and fl m is the density parameter 
of non-relativistic matter (e.g., baryons and/or weakly interacting massive particles). We 
will focus on recent epochs, when the radiation energy density is negligible compared to 
non-relativistic matter. 

The cosmic evolution of <fi is essentially determined by the ratio of its mass, ~ M 2 //, 
to the instantaneous expansion rate, H(t). For ^ 3H, the field evolution is overdamped 
by the expansion, and the field is effectively frozen to its initial value. Since <f> is initially 
laid down in the early universe (at a temperature T ~ / 3> M) when its potential was 
dynamically irrelevant, its initial value in a given Hubble volume will generally be displaced 

7 



from its vacuum expectation value (f> m = irf (vacuum misalignment). Thus, at early times, 
the field acts as an effective cosmological constant, with vacuum energy density and pressure 
ptf, ~ — pj, ~ M 4 . At late times, m$ 2>> 3H(t), the field undergoes damped oscillations 
about the potential minimum; at sufficiently late times, these oscillations are approximately 
harmonic, and the stress-energy tensor of <f> averaged over an oscillation period is that of 
non-relativistic matter, with energy density p$ ~ a~ 3 and pressure p^ ~ 0. 

Let t x denote the epoch when the field becomes dynamical, m$ = 3H(t x ), with cor- 
responding redshift 1 + z x = (a(ta)/a(t x )) = (M 2 /3i?o/) 2//3 ; for comparison, the universe 
makes the transition from radiation- to matter-domination at I0 4 £l m h 2 [where 

h = i?o/(100 km/sec/Mpc)]. The f — M parameter space is shown in Fig. 1. To the right 
of the diagonal line = 3i?o, the field becomes dynamical before the present epoch and 
currently redshifts like non-relativistic matter; to the left of this line, <f> is still frozen and 
currently acts like a cosmological constant (the region denoted by 'A'). In the dynamical re- 
gion, the present density parameter for the scalar field is approximately ~ 247r(//Mp/) 2 , 
independent of M [12] (assuming the initial field value = 0(1)/); thus, the horizontal line 
at / = 1.4 x 10 18 GeV indicates the cosmic density limit $1^ = 1. In the frozen (A) region, 
on the other hand, fl^ is determined by M 4 , independent of /, and the bound Cl^ = 1 is 
indicated by the vertical line. 

Focus on the dynamical region in the right-hand portion of Fig. 1. If <j) dominates the 
energy density of the Universe, the growth of density perturbations is strongly suppressed 
for physical wavenumbers larger than the 'Jeans scale' [16] kj ~ m^(<^ m (£) / Mp/) 1//2 , where 
$m{t) ~ /[(l + - 2 (^))/(l + z x)] 3 ^ 2 is the amplitude of the homogeneous field oscillations 
at z(t) < z eq . If this Jeans scale is too large, perturbations on galaxy and cluster scales 
would not grow at high redshift, leading to a power spectrum with an unacceptably large 
coherence scale. We can express the resulting perturbation power spectrum in terms of the 
standard cold dark matter (CDM) spectrum as P(k) = P c dm(k)F 2 (k)] for z x > z eq , the 
relative suppression factor due to the scalar field is [17] 
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Fig. 1: The PNGB model parameter space. 
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Here, 1 + z*(k) = [(M/A;)(3i?o/Mp/) 1/ ' 2 ] 4 is the redshift at which the physical wavenumber 
kphys = k(l + z ) drops below kj, so that scalar perturbations on that scale can begin to grow. 
Thus, M sets the scale where the power spectrum turns down from the CDM spectrum, and 
/ (through ft,!,) determines the spectral slope n of the suppression factor, F(k) ~ k~ n with 
—4 < n < (note that for fl^ ^ 0.2, n ~ 1241^/5) . For galaxies and quasars to form at 
moderate redshift, the power at small scales should not be very strongly suppressed compared 
to standard CDM. We therefore impose the appoximate bound F(k = 1.6/iMpc _1 ) > 0.3, 
which corresponds to the curved boundary in Fig. 1: the region above this curve is excluded. 
To the right of this region (in the area marked CDM), <f> acts as an ordinary cold dark matter 
candidate, a lighter version of the dark matter axion. In the area marked MDM, the effects of 
<j) on the small-scale power spectrum are similar to those of a light neutrino in the mixed dark 



9 



matter model: at the point marked by the star, the variance of the density field smoothed 
with a top-hat window of radius R = Mpc is <7g((j)) ~ <jg dm /2. When the amplitude is 
normalized to COBE on large scales, this yields crg(0) ~ 0.6, as suggested by the abundance 
of rich clusters of galaxies and the small-scale pairwise velocity dispersion of galaxies. In 
this region of parameter space, the neutrinos of mass m v ~ M ~ several eV could play a 
dynamical role in structure formation as well. 

For the remainder of this Letter, we focus on the parameter region near the bullet in 
Fig. 1, in which the field becomes dynamical at recent epochs, z x ~ — 3, or in the near 
future: this has new consequences for the classical cosmological tests and the expansion 
age, and it does not lead to the small-scale power suppression above. We thus impose the 
constraint m$ = M 2 / f ^ 3Hq. The second condition is that the PNGB energy density be 
dynamically relevant for the recent expansion of the universe, which implies p^(^o) ~ pcrit(ta), 
or M 4 ~ 2>H 2 M 2 n l&Ti. Combining these two constraints determines the two mass scales 
in the theory to be / £ M p/ /(24tt) 1 / 2 ~ 10 18 GeV and M ~ 3 x lO" 3 /* 1 / 2 eV. As argued 
above, we can construct particle physics models for light PNGBs with these mass scales: the 
spontaneous breaking scale / is comparable to the Planck scale, and the explicit breaking 
scale M is comparable to that expected for light neutrinos for the MSW solution to the 
solar neutrino problem. The mass of the resulting PNGB field is miniscule, ^4x 10~ 33 
eV, and (by construction) its Compton wavelength is of order the current Hubble radius, 

= m^ 1 = Hq 1 /3 <: 1000/i -1 Mpc [18]. This is a generic feature of scalar field models for 
relic vacuum energy that satisfy V((j) m ) = 0. 

Figure 2 shows several examples of the evolution of the scalar field [Eqn.(7) with the 
potential of Eqn.(6) and the Hubble parameter given by the expression immediately below 
Eqn.(7)]. We show £l m = 1 — £1^ as a function of the expansion age Ht, for different initial 
values of the field <f>i/ f (assuming = 0, since the field is Hubble-damped at early times). 
The numerical evolution starts at p m /M 4 1, i.e., at the top of the figure (Q m ~ 1 Q^) 
in the matter-dominated epoch. At early times, the field is effectively frozen to its initial 
value by the Hubble damping term in Eqn.(7), and the evolution tracks that of a cosmological 
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constant model (curve labelled 'vac' in Fig. 2). At t ~ t x , the field begins to roll classically; 
on a timescale initially comparable to the expansion time, the expansion age Ht reaches a 
maximum and subsequently falls toward 2/3 (indicated with the vertical dashed line) as the 
field undergoes Hubble-damped oscillations about the potential minimum. The evolutionary 
tracks are universal: a shift in the mass scale / accompanied by an appropriate rescaling 
of the initial field value & leads to essentially identical tracks, i.e., a given track actually 
corresponds to a family of choices of (^, /). 




0.6 0.8 1 1.2 

Ht 

Fig. 2: The non-relativistic mass density fl m = 1 — fl^ vs. Ht, for / = Mpij 'V87T '. The solid curves 
correspond to several initial values for the field, = 1.4, 1.5, 1.6, and 1.75. The evolution starts at the 

top of the figure and ends at the lower left. The vertical dashed line shows the Einstein-de Sitter expansion 
age Ht = 2/3, the horizontal dashed line shows the lower bound Q m = 0.1 from dynamical mass estimates, 
the dotted curve (labelled 'vac') shows the evolution for a cosmological constant model, and the long-dashed 
curve corresponds to an open model with Q,^ = 0. The dot-dashed curves (labelled 0.6, 0.7, 0.8) bracket the 
constraints from lensed QSOs in the HST snapshot survey (see text). 

The observational consequences of this model follow when one identifies the present epoch 



11 



to on an evolutionary track — this implicitly corresponds to fixing the mass scale M. For a 
given expansion age HqIq, one can choose the upper branch, where the field is still frozen and 
thus nearly identical to a cosmological constant, or the lower (dynamical) branch, for which 
the recent evolution will be intermediate between vacuum- and matter-dominated and which 
has qualitatively new features. Dynamical estimates of the mass in galaxy clusters indicate 
the lower bound fl m j£ 0.1 for the mass density in non-relativistic matter. Consequently, the 
lower branch is excluded if the initial value of the field is below some value, e.g., <f>i/ f — 1.3 
for / = Mp//\/87r. Physically, for such small values of <f>i/ f, the universe undergoes several e- 
foldings of inflation before the field begins to oscillate, diluting the density of non-relativistic 
matter. Consequently, to achieve large expansion times in this model, HqIq ~ 1, the present 
epoch must be in the vicinity of the 'nose' of the evolutionary track, which corresponds 
approximately to the condition t x ~ to imposed above. 

As with vacuum-dominated models, these scalar field models can in principle reach ar- 
bitrarily long expansion ages, Ht ^> 1, if 4>il f ls sufficiently small. However, this region 
of parameter space is excluded by the observed statistics of gravitationally lensed quasars. 
The 3 dot-dashed curves in Fig. 2 show the observed constraints on the incidence of lensed 
QSOs. We computed the number of lensed QSOs expected in the HST Snapshot survey [19] 
for cosmological constant models with fl\ = 0.6, 0.7, and 0.8; along the 3 curves in Fig. 
2, the number of expected lensed QSOs in the PNGB models are equal to these 3 values. 
Since different assumptions about galaxy models yield different lensing fractions, we show 
the limits corresponding to these three cases to cover the spread of quoted limits in the 
literature [7] (the region to the right of each curve is excluded). For a given lensing limit, 
the upper bound on the expansion age HqIq is increased in the scalar field models compared 
to the cosmological constant model; imposing the lower bound £l m > 0.1, the bound on i?o^o 
can be relaxed by 7 — 10%. Thus, the scalar field models are relatively more successful than 
a cosmological constant at easing the 'age crisis' while remaining within the observational 
constraints, provided fl m is fairly low. 

We have presented a class of models which give rise to (technically natural) ultra-light 

12 



pseudo-Nambu-Goldstone bosons. With spontaneous and explicit symmetry breaking scales 
comparable to those plausibly expected in particle physics models, the resulting PNGB 
becomes dynamical at recent epochs and currently dominates the energy density of the 
universe. Such a field acts as a form of smoothly distributed dark matter, with a stress 
tensor at the current epoch intermediate between that of the vacuum and non-relativistic 
matter. Such a model 'explains' the coincidence between matter and vacuum energy density 
in terms of particle physics mass scales, reconciles low dynamical mass estimates of the 
density parameter, fl m ~ 0.2, with a spatially fiat universe, and does somewhat better than 
a cosmological constant at alleviating the 'age crisis' for spatially flat cosmologies while 
remaining within the observational bounds imposed by gravitational lens statistics. 

This work was supported by the DOE and NASA grant NAG5-2788 at Fermilab. JF 
thanks Lloyd Knox for his scalar field evolution code and Richard Watkins for discussions. 
After this work was completed, we became aware of related work by Fukugita and Yanagida 
[20], which considers an axion model for the non-dynamical (A) region of Fig. 1. 
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